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Abstract: We study preheating in DBI hybrid inflation. Preheating occurs since the 
matter fields are non-minimally coupled to the inflaton. Despite the coupling being small, 
preheating happens as a narrow parametric resonance in which matter fields and inflaton 
modes are exponentially created. We show that, within the approximations used, the 
process is very efficient at late times. The resonance stage gives way to the decay of 
the inflaton into bosonic or fermionic degrees of freedom and the corresponding reheating 
temperature is determined. 
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1. Introduction 

Recent developments in string theory have resulted in attempts to construct string inspired 
models of inflation. Whilst progress has been made in brane inflation models ||l|, ^, ^, the 
theory of reheating at the end of inflation in such models is still in its infancy 0, |^ • 
Most string inspired models of inflation have concentrated on the form of the potential. 



However, in an interesting paper [10| Silverstein and Tong proposed a novel D-cceleration 
setup within the context of Brane Inflation theory. This opened the possibility for fast- 
roll inflation to take place, circumventing the fine-tuning problems. Their model uses 
non-standard kinetic terms for the inflaton field of the Dirac-Born-Infeld {aka DEI) type, 
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motivated since the DBI action provides the low-energy description of a D3-brane in an 
AdS space |11]. 

A novel feature of DBI inflation is that it predicts non-Gaussian primordial spectrum 
12] for the scalar perturbations which can then be confronted with observational data. 



In addition, the speed of sound is time dependent [13|. Available data suggests tighter 
constraints on the parameters of the brane inspired inflationary models. This opens up 
the exciting possibility for the Planck satellite |14] [^] to shed some light on other models 



besides the standard single-field inflationary ones [^]. Whereas the inflationary perturba- 
tions have been studied in DBI inflation, up until now there has been no attempt to study 
the theory of reheating. In this paper we study a simple model of reheating within DBI 
inflation. 

Reheating occurs at the end of inflation when the energy stored in the inflaton field 
is transferred to matter and radiation due to the oscillation of the inflaton around the 
minimum of the potential. Due to the coupling with matter fields this oscillation can lead 
to parametric resonance, resulting in the explosive non-perturbative growth of particles. 
When the final state is thermalized, it is said that reheating has occurred. This preheating 
phase has been studied extensively for standard models of inflation ||l^, ^ |6|, |l8|. The 
extension to models with non-standard kinetic terms is in its early stages, though non- 
standard matter fields have been studied recently |1£]. In this paper we take a different 
approach. We consider DBI inflation with the inflaton field coupled to matter fields. As a 
first step we take the matter fields to have standard kinetic terms. 

Our model is a variant of hybrid inflation ||2^ and is motivated by D3/D7 brane 
inflation. We show the modifications to the preheating effect that need to be taken into 
account, when the inflaton field has DBI action. The theory is described by 

5* = y V-9m ' (Co + ...) , 
with signature (— , -|-, -|-, -|-) for the metric and 



Co 



V{cP,x)-i^d^xd''x, (1.1) 



where the inflaton field (p is measured in units of the string coupling, gg. f{(j)) is the warp 
factor, which is inversely proportional to the brane tension and x is the matter field. The 
potential used is 

V ^V{^,x) = 2gW + '^{x'-C3)' , (1.2) 

where ^3 is the symmetry breaking term. Any one-loop corrections are sub-dominant and 
will not be considered. 

The modifications due to the DBI kinetic term will be studied and compared to the 



case with standard kinetic terms and the same potential |21|. This way we gain insight 
into how (p)reheating is modified due to non-standard DBI kinetic terms. This paper is 
organised as follows: first we study the cosmological evolution of the inflaton field as it rolls 
down the potential from the false to the true vacuum configuration. We then assess the 
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differences in tfie (p)relieating process, analysing the possibility of broad/narrow resonance 
occurring. We compare the efficiency with respect to the usual case. Finally, we discuss 
our results and we summarise the underlying implications. 

In what follows we will use natural Planckian units, in which the reduced Planck mass 

is npi := ~ 2.4 X lO^^GeV. 

vSvr 



2. Dynamics of Inflation 

In this section, we compute the cosmological evolution of the scalar fields and investigate 
the circumstances under which (p)reheating occurs. 

2.1 Equations of Motion 

The equations of motion for both fields are 

* + 3/f;^-i5 + -^.*^ + ;^(I, + ^j =0, (2.1) 

and 



Here 



dv 

X + 3i/x + = • (2-2) 



7:= (l-/</>^) (2.3) 



plays a similar role to that of the Lorentz contraction factor and dot denotes the derivative 
with respect to cosmological time. Current data requires 7 < 32 for compatibility with 
observation p^ . 

Clearly, the difficulty with these equations is that the square root from the inflaton 
kinetic term makes the equations highly non-linear. To solve the equations we adopt 
approximation techniques. Firstly, we assume a cut-off throat in the infrared, so that the 
warp factor becomes 

/(</,) = , 

for some new energy scale [i. Since we are interested in the late time behaviour towards 
the very end of inflation, then (/> ^ 0, which allows us to simplify the warp factor as 

m) = ^ . 

This assumption will be in agreement with the further approximations in the next sections. 

The modiflcation of the equations of motion for the inflaton field result in a modifica- 
tion to the the usual preheating mechanism. In particular, the solutions to the equations 
of motion of the inflaton field acquire a correction to the canonical harmonic oscillator 
equation. We make this explicit in the next section, where such (perturbative) correction 
is computed. 
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2.2 Reaching the bottom of the potential 

At the start of inflation, (x) = and the inflaton starts roUing down the potential. As 4> 
goes down the potential, slow-roll conditions are no longer valid (since the field speeds up 
when meeting a steeper region of the potential) and it approaches its bottom, overshooting 
and it starts oscillating. Eventually, both fields should assume their vacuum expectation 
values, namely 

(</.)= and (x) = V^ , 

for which they minimise the potential and violate slow-roll conditions. 

Next, we will study the fluctuations of the flelds around their background solutions. 



3. Preheating 



In our study we neglect the Hubble friction term since the expansion time is much larger 
than the one we are interested in. We perform a Taylor expansion of the 7 factor, keeping 
first order terms, and thus neglect O { f"^ 4>^] since this is higher order. This is because 



the presence of the branch cut in the denominator of (2^) demands a limiting speed of the 
D3-brane moving along the throat to be 



< 



(3.1) 



for (p ^ The modulus sign accounts for the fact that (p decreases in time as it gets closer 
to the tip of the throat. Hence, the contributions coming from these higher order terms 
are negligible. Then, ( |2.lD becomes manifestly simplified and 



+ 



'-4 

2 







(3.2) 



Considering the extra term to be a small perturbation to the homogeneous harmonic oscil- 
lator equation, we require f cfP' < 2/3, which is a stronger requirement but still compatible 
with (p.l|). We hence get a correction to the usual preheating equation of motion for the 
inflaton, with solution given by 

(f) = A cos {oj t) , 

with A and uj = 2g^/^ being the oscillation amplitude and frequency, respectively. It 
follows that the solution to ( |3.2D may then be obtained by assuming an additional pertur- 
bative contribution to the latter equation, yielding 



3 A 

^cos {ujt) + - —A^oj'^ sin^ (ujt) cos (ujt) 



(3.3) 



up to a normalisation factor. This ansatz is still a periodic function of time, so that one can 
trace its time evolution throughout an unlimited number of periods. The strength of this 
perturbation is parametrised by the coefficient sitting outside the periodic function and we 
note that, due to the Hubble friction term, the amplitude of the oscillations is expected to 
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3 A 

be damped. Thus, we expect — jA^lj'^ < A, which is in agreement with the condition for 

4 /X* 

the stronger restriction < 2/3. 

Slow-roll conditions stop being valid when, for small values of the matter field X) 



< 



— and a tachyonic instability sets in. From here and to lowest order in the 



approximation scheme, the oscillatory background solution for the matter field 
follow a linear trajectory according to 

X=VC3+X with x = -V2(t) , 



will 



(3.4) 

which can be shown to obey the equations of motion and ( |2.2| ). This suggests the 

correct rescaling in ( |3.3D for the inflaton background solution. Thus, we consider the 
(normalised) classical background solutions to the equations of motion ( |2.1D and ( |2.2D to 
be 



A 



\/2'(^ = ^cos (tot) + 3g'^ -T-^^Cs sin^ i^t) cos {tot) 



(3.5) 



X 



^ - A cos (ujt) - 3 5^ ^-^^Cs sin^ {ujt) cos (wt) . 



3.1 Fluctuations around the background solutions 

We now consider the approximate classical fluctuations in the direction of the background 
solutions in ( p.Sp , defined by 

6(f) := (f> — (f) 

(3.6) 

^Sx-=x-x ■ 

Importantly, in order to solve the modified equations of motion ( ^.l]) and (^]^) we need 
to consider O {A^) ; however, for the purpose of studying the equations of motion for the 
field fluctuations it is sufficient to consider O (^^) terms. Also, we consider the fluctuating 
Fourier modes in a Hubble volume with comoving number k (up to a normalisation factor): 



) = d^k {ar(j)kit)e 



Xix) = d^k {brxkit)e 



-ik.x 



+ aicf>l{t)e 



ik.x 



+ blxUt)e 



ik.x 



where {ag, 6^} and I'^^j^gl respectively the annihilation and creation operators, with 

momentum k. Moreover, since we are expanding around the background solutions (p.5|), 
which satisfy the equations of motion (2.1) and (p.2|), it follows that only the quadratic 
terms in the potential ( |1.2D are relevant. Hence, the fluctuations with comoving number k 
obey the following equations of motion in the momentum space: 

3 A 



+ k^dcpk + 



1-- — </,2 



85^ 4>X^Xk+ V X^^k 







(3.7) 



Sxk + Sg"^ (j)x5(t>k + 



k' + 2g' ( 2^' + 3x^ - Cs 



6xk ^ 
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In particular, in agreement with the hnear trajectory defined in ( ^.41) , the fluctuations of 
both fields are related by 

bXk = -^^4>k . (3.8) 



Plugging it into (|3.7D and using the ansatz (^) , we get the simplified equations of motion 
for the fluctuations of the inflaton and matter fields around their respective background 
solutions, up to the order of approximation used: 



h + 



A 



4 y ''3 



+ 



-Wg^ \/C3A cos (tut) + 6g^ [1 + ^g^Ci] cos (2 ujt) 



(^Xfe + [fc^ + VCs + 65^2 - 16^2 cos (cjt) + Qg'^A'^ cos {2ujt)] 5xk - . 

(3.9) 

We note that the equation for the fluctuations of the inflaton field resembles Hill's equation, 
which is usually written in the form 



00 + 2^ COS {2rz) 



r=l 



y = o , 



+ 00 



where 6*0, O2, 64, ... are assigned parameters and |02r| converges. Indeed, introducing 



r=l 



a rescaled dimensionless time defined as 2z := cut, it follows that the first equation in (| 
may be written as 



+ [60 + 26^2 cos (2z) + 26*4 cos (iz) 







(3.10) 



where the prime denotes the derivative with respect to the rescaled time with the identifi- 
cations 



A 



< 



A 



(3.11) 



Hill's equation with three terms may be transformed in a type of equation which was 
already studied in [19| by performing the following change of variables |23| 

5(^k = Vke~^ cos{2z) , 

so that in terms of the new field, ( 3.10| ) is rewritten as 



^ + 4V^sin (2z) ^ + 60 + 264 + cos (2z) (262 + ^^/H 



-6- 



This has the generic form of 



^ + 2psm(22)— + [i4 -2^cos(2z)]^; = , 



(3.12) 



with 



V 



+ 29^ 



(3.13) 



KQ =-{92 + 2^07). 
We observe that q can be either positive or negative, depending on the relative strength of 



62 and 64. This equation is known in hterature |24] and we will study this in detail in the 
next section. 

It is intriguing that the phenomenon of inflaton particle production also takes place 
in preheating. This follow as a consequence of the DBI term. In fact, at the level of the 
equations of motion, ( p. 21) may be rewritten in a more suggestive way as 

1 i ^^^// 



with Veff approximately given by 



eff 







(3.14) 



(3.15) 



and with V as defined in (|1.2| ). Therefore, the rescaled potential has now a cubic power in 
the field, explaining the self interaction observed here as a consequence of the DBI term. 
In this sense, DBI models fall in the same class as self-interaction models where resonant 
infiaton production is, in general, expected [p^ j. 

We note that since the corrections to the inflaton fluctuations are we can 

truncate the contributions to both the matter and inflaton equations at this order. This 
amounts to neglecting the DBI correction terms to the inflaton background solution. There- 
fore, we would expect to obtain the same equations of motion for the fluctuations about 
the background matter fleld solution as in |21|, except that now O (-4.^) can no longer be 
neglected as being subdominant, a priori, at this order. However, the matter field equa- 
tions involve no corrections from DBI^. This allows us to ignore these terms, obtaining the 
further approximate equation of motion: 



16(7 V Cs cos (wt) 



5xk ^ 



(3.16) 



This has the standard form of Mathieu's equation which is a special case of Hill's equations 
by taking ^4 = in ( 3.11 ). Indeed, introducing the dimensionless time scale defined via 
2z = ojt and denoting, like before, the derivative with respect to 2; as a prime, then the 
above equation is rewritten as 



5xk + [Ak-2q cos {2z)] <5xfc = 



(3.17) 



^this is because in the case of the fluctuations of the inflaton fleld the strength of the terms is not only 
parametrised by A but also by the cut-off radius ^ in the term A//^*. 
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where 

g%3 

(3.18) 

Mathieu's equation is well known (see, for example, |^ and ||2^) to be associated with 
preheating. Its features are essentially that its solutions are a linear combination of expo- 
nentially growing and decaying modes, leading to instability bands dependent on the choice 
of values for and q, which are in turn interpreted as a particle production phenomenon. 

3.2 Efficiency of particle production 

Although the equation for the matter field fluctuations in (|3.9|) was already obtained in 



1 21], the equation for the inflaton field fluctuations contains new terms, requiring a concise 



study of the orders of magnitude of the terms in ( 3.10 ). The main point is that the speed 
limit (|3.1| ) for the D3-brane moving along the cut-off throat provides a natural limit for 
the orders of magnitude of several parameters in this equation. Together with the lowest 



order background solution for the inflaton field (|3.5|) , it follows that 

/x^ 3 

This, together with the decreasing amplitude of the oscillatory contributions themselves, 
verifies that the correction to the usual background solution in ( |3.3| ) is of lower order. In 
order to obtain an order of magnitude we can take the average of the oscillatory solution 
over several periods with respect to the late time behaviour. Since y/Cs' and fi are fixed 
mass scales, the order of magnitude of the term in the criteria will vary as A decreases 
in time. Therefore, its maximum value corresponds to the maximal oscillation amplitude 
A in this regime. However, at early times we cannot use the same criteria for averaging 
over periods. However, if /^^ < 2/3, then this will certainly be true when the scalar field 
has maximal speed, ie, at the very bottom of the potential. Furthermore, by the time the 
field has reached the bottom of the potential during the oscillation, the amplitude A has 
not changed significantly. Taking this into account, the following criteria should follow, for 
early and late times, respectively: 

For now, "late" times means studying the oscillation over enough periods so that the 
averaging over the oscillatory solution is approximately (sin^ (i^^)) — 1/2. 

To proceed, we need to compare the two energy scales, -v/Cs" and A. Following for 
the usual kinetic terms, the condition for efficient resonance was shown to require that the 
final amplitude Af for the oscillations of the field around the minimum of the potential 
were related to the initial amplitude Ai by 

^/ < f^V" . (3.20) 



Ai V64/Upi 
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Af 1 

To get an order of magnitude, we assume the limiting value for which —r- ~ — . This will 
be tested later against the condition for efficient preheating. Hence, at the very beginning 



of the oscillation, since we expect Ai ~ vCs > then 



A^ 



1 



on the other hand, as a consequence of particle creation taking into account the contribution 
of the quartic self-interaction term in x in draining energy from the inflaton field, we expect 
A to decrease significantly, so that 



A^ 



1600 



A 



Thus, it follows that the term —g C3 becomes relevant at late times, when the DBI 

corrections are more important. 

In this spirit, we may approximate the equations of motion for the inflaton modes 
in (|3.9|) into two different regimes: up to a few oscillations and then the truly late time 
behaviour. 



3.2.1 Within a few oscillations 

At early times, the inflaton mode equation in ( |3.9D may be approximated by 
5<t>k + Ik"^ + 45^(3 + %^ A^ - 165^ v^-^ cos (wt) + 65^ A"^ cos (2 ujt)] 6(pk ^ 



, (3.21) 



for Af'^A< VCs- 

In this regime we can make further approximations by taking into consideration which 



terms contribute the most to this equation. In particular, because 



equality will only be made manifest at the beginning of the oscillation. Hence the effect of 



the quadratic terms in A may be neglected so that for A < \/0r/4, (|3.21 ) takes the form 



k + 



I65 V Cs cos (ujt) 



~ 



(3.22) 



We recover the same equation as in |21|, so we only state their main results. This is also 
equivalent to ( 3.17 ), so we will be seeking unstable solutions of this Mathieu's equation 
(this is as expected as the DBI corrections become important at later times). This solution 
is of narrow-resonance type, but is efficient since the modes remain in the resonance band 
for sufficient time to become exponentially modified. In fact, denoting the instability 
parameter with no DBI corrections by the result is that the unstable mode should 
grow as 

(3.23) 



JgAt 



we identify the width of momenta of produced particles as A/c = 16/3\/3'gA with centre 
at ko = g\/3 Cs 
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3.2.2 At later times 

After a few oscillations, we go into the late times regime for which (p. 91) takes the form 



S4>k + (k" + 4^2(3 - Q\g^A\l + Q\g^A^Cl cos (2u;t) ]5<Pk^^ 
\ II ji 

which is of the form of Mathieu's equation 



+ [Ak- 2^cos(2z) = , 



where 



Z = LOt 



(3.24) 



(3.25) 



1 + 



_ 3 ^ 



(3.26) 



The criteria for narrow resonance is verified [26| 

» 9'C3-l^,9'C!A' > l^.9'CiA' . 

This essentially corresponds to Ak > 2 The interesting behaviom' of Mathieu's equation 
is due to the presence of two master frequencies: Ak and |2gcos(2z)|. The interplay of 
these two frequencies, particularly since the second one is evolving in time, generates 
interactions between the two possible regimes depending on whether A^ < |2gcos(2z)| or 
Ak > |2gcos (2z)|. 

The theory of narrow parametric resonance is well known ([^ and In short, the 
resonance occurs for very narrow bands in the momentum and the most important and 
widest is the first one. Indeed, for such case, Ak ~ 1 ± |q| and so the values of k which are 
affected by the resonance stage will be lying within the range 

The general solution to Mathieu's equation may be written as 

6M^) ~ ae^'^'fiz) + be-'''^'f{-z) , 

where ^k is the instability parameter (slowly varying function of z) and f{z) is some periodic 
function in z with period tt. Solutions with 3^(^) ^ are superposition of exponentially 
decaying and growing modes and these last ones may induce resonant behaviour. For the 



narrow regime, using Hill's recurrence perturbative method |23], one can show that the 



maximal value for fj,k is Hkl 



\q\ /2, ie, 



(3.27) 
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This implies that the growing modes behave as 

6(j)k ~ e**^ , which maximised yields Scj)^ ^ e 2 ; 

the maximum value will happen for the maximal speed of the brane moving along the 
throat. This results in particle production. 

Comparing the two maximised instability parameters from the early (3.23) and late 
( 3.27 ) time stages, we verify that fik > f^k/"^ some oscillation amplitude, A, of the 
inflaton field. This happens beyond a certain time. We conclude that the DBI correction 
boosts the explosive particle production when 

A<A=^^ = 0.1257C^ , 

ie, before the end of the oscillations, when A = Af ~ 0.025v^^. 

This behaviour should be true for adiabatic changes in /x; otherwise, the growing mode 
would be of the form 

However, we can compare the rate of the time evolution of ^ with the typical frequency in 
our theory, w, as the suitable criteria for determining the adiabaticity of the resonance, 

— < w ; 

this is obeyed when the time evolution for the oscillation amplitude. A, is such that 

A^A{t)^^ . 

Provided this is true, then the resonant modes behave as described above. We expect A to 
be decreasing in time not only as a consequence of the Hubble friction term, but also due 
to the quartic interactions of x which contribute to both rescattering and back-reaction. 
These factors all contribute to reduce the resonance effect. 

3.3 Is preheating efficiency undermined by the narrow resonance type? 

Since the resonance is narrow, rather than broad, we need to test whether or not it is 
efficient. Thus we need to compare the time each mode stays inside the resonance band. 
The necessary criteria is 

f2uj>H . 

This condition is verified for all amplitudes A <^ Ac 0.Q5 fipi ^ 220 ^/ (3, which is the 
case for expected values of -v/CF 1 21]- Instead, if we used the Hubble parameter at the end 
of a slow-roll period 

rr 9 CS 
11 — 

then the criteria above would be satisfied for \/C3 ^ fJ-pi, in agreement with the 
standard case. 
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4. Reheating 



Due to the explosive behaviour of the resonance stage, the matter particles created via 
this process are far out of equilibrium. As previously discussed, not only matter particles 
but also inflaton particles are produced via the parametric resonance stage. Since this 
can have dangerous consequences, we need to provide a means for the inflaton particles 
to decay naturally into matter particles. Thus we investigate possible decay channels to 
obtain standard reheating. In this section, we also comment on the energy scales involved 
and how this affects the reheating temperature. 

4.1 Power spectrum constraints 

To proceed we first constrain the parameters of the theory, namely the symmetry breaking 



and cut-off energy scales normalised with respect to the CMBR |27, For standard 
kinetic terms, the density perturbations are given by 

- - • ("> 

P 7r3/2 (j) 

for a time evolving classical scalar field. For DBI models, due to the normalisation of the 
scalar field with respect to the string coupling, this becomes 

p 7r3/2 (f) 

Since density perturbations are generated at the end of inflation and before preheating 
only the inflaton fleld will be considered. In the regime for which f cfp' < 2/3, the maximal 
speed of the inflaton fleld is 

1 PT 






< 











\/rv3 • 

Since the inflaton field is rolling down the potential, we take the negative root. At this 
stage the potential is 

plugging in Friedmann's equation, we get 

= ' = 1 «M . (4,3) 



Thus, matching with COBE data, it follows that 

:V)' = ^A4.4xlO-'» . (4,4) 

p ) 247r'^ gs p^ 

In it was estimated that if the inflaton is to drive a purely de Sitter inflation, then 
this would lead to bounds on the 't Hooft coupling as A/gf<j ~ 10^^. Thus we expect, for 
reasonable values of 5, Cs ^ (a*) Ppi'i this is consistent with our choice for the warp factor. 
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4.2 Available decay channels 

Our aim is to propose a scenario in which decay channels are available as soon as the 
resonance stage is over so that the resonance process can be shut down at late times. On 
the other hand, these should not undermine it, which immediately discards the possibility 
of the inflaton decay via (jxj) XX- Ultimately, as a natural consequence, we consider that 
any massive scalar matter particles produced here will eventually cascade down into the 
Standard Model particles and hence reproduce the content of our Universe. 

The proposed channels through which the residual energy density of the inflaton is 
transferred to matter are only efficient if 

(x^) > i'P^) 1 ^ i'P) ) for the cubic channel 

^ ' ^ ' 2 zg-^ 

and 

{x^) > {4>) , for the fermionic channel, 
which is true when the fields reach their true vacuum values. 
4.2.1 Cubic interaction 

We consider adding new interactions to the Lagrangian density Cq of the form X^cpx^, 
where is the coupling constant, with dimensions of mass. It can be shown that the 
introduction of this term does not affect the preheating analysis provided 

^ — « voT- 

2 

We shall consider such regime. For this process, the decay rate is given by B 



r(0^xx) = ^. (4.5) 



This should be compared with the resonance rate 



^ ^ _ ^ ^o; . (4.6) 



Using ( 3.271 ) we evaluate the condition T > fiuj for which the amplitude is 



gs 127r A g%3 

For the final amplitude Af = \/C3'/40 (in gg units), the resonance turns off, and using 
( |4.4| ), we get the minimal value for the coupling constant for this decay channel: 

A^ > (7 X 10-9) 4^ . (4.8) 

Both the estimates in ( |4.4D and (^^) will be relevant when evaluating the reheating tem- 
perature. 
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4.2.2 Adding fermions 

We could however expect other degrees of freedom to be present in additional decay chan- 
nels. Thus we consider the coupling of the inflaton field to fermions. Although this is 
less efficient as far as the preheating mechanism goes, we can still evaluate the reheating 
temperature. We add the term —hipipcj) to the Lagrangian density, where '0 denotes the 
fermion field and /i is a dimensionless coupling constant. The decay rate is given by 

r(^^V^V) = ^ • (4-9) 

If this is to prevail over the resonance stage at late times, then 

A"^ h? 1 
gs Svr A 5^ Ca 

4.3 Reheating temperature 

We can calculate the reheat temperature from the energy density. The energy density in 
particles is 

where Te// is the effective decay rate. However, for relativistic particles in thermal equi- 
librium 

p^'^n{Tr)tI, , 

where the number of degrees of freedom, N{Tji), is of order 10^ or 10^ at the reheating 
temperature, Tr. Thus 

Tr ^ 0.2 ^FeffMpi' . (4.11) 

4.3.1 Inflaton-Matter interaction 

Plugging in ( [4. 5]) , then 

Tr^ 2.1xlO-2Af^y ' . (4.12) 



From the COBE bound (|4.4| ), we get an expression for the symmetry breaking scale 
as a function of the coupling g and the cut-off scale ^: 

(3=^1.65x10-5^ . (4.13) 
9 

Using ([4.8| ), this gives the reheating temperature to be 

5/2 

T, :.5xl0--^ . (4.14) 

Inserting the constraint for the symmetry breaking scale in ( 4.13| ), we finally get an estimate 
for the reheating temperature by eliminating the (3 dependence: 

7/4 

Tr ^ 2.1 X 10- ^ . (4.15) 
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So, for g = 8 X 10^ , then a reheating temperature of T ~ 7 x W °GeV corresponds to a 
cut-off energy scale of ~ 9 x 10~^/ip;. A general feature of ( ^.15| ) is that for more warped 
throats the reheating temperature will be considerably larger. We note that for standard 
kinetic terms and coupling constant g ~ 10~^, T/j ~ 4.2 x lO^^GeV. Thus, DBI inflation 
seems to predict a higher reheating temperature. 

4.3.2 Infiaton-Fermions interaction 

Proceeding in a similar fashion for the fermionic decay channel, taking X/gs ~ 10^^, we 
find that the reheating temperature is 

r« -29.1 5^/2 (^^j ^.P, . (4.16) 

For g ~ 10~^, one obtains Tr ~ 6 x lO^^GeV for /z ~ 9 x 10^^ Mpi, which is again quite 
high. 

5. Conclusions 

We have studied preheating and reheating in DBI inflation, coupling the inflaton to matter 
fields and using a hybrid potential. The expansion of the Universe was neglected in this 
work. By considering the DBI terms to be a small correction to the standard kinetic terms, 
we have determined the cosmological evolution of the inflaton as it oscillates around the 
minimum of the potential; we have also studied the evolution of the matter field. We 
have shown that preheating proceeds in the narrow resonance regime, which proves to be 
efficient in draining energy out of the inflaton field. Interestingly, DBI correction terms 
enhance the instability associated with the comoving modes of both fields after several 
oscillations. This suggests that non-standard kinetic terms favour efficient preheating by 
allowing stronger resonance of both the inflaton and the matter modes. 

In this paper, it was proposed that preheating was followed by the elementary theory 
of reheating. The conditions for this to happen constrained the coupling between the fields. 
Current bounds from the CMB spectrum were used to obtain estimates of the reheating 
temperature. This is a crucial step in order to achieve a perception of the energy scales 
in the theory. General estimates suggest that DBI inflation models give larger reheating 
temperature compared to theories with standard kinetic terms. 

Although the work presented here is based on a toy model, it should give some insight 
into preheating in DBI inflation. We regard this as a first step. Future work should consider 
the expansion of the Universe and go beyond the approximation used in this paper. It 
would be interesting to draw a general conclusion of DBI inflation in (p)reheating with 
more general potentials [pS|]. 
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A. Classical Dynamics considerations 

We would like to obtain the asymptotic behaviour that (/> — > 0, at the bottom of the 
potential, using Classical dynamics principles. We intend to get such behaviour before (j) 
reaches the true vacuum and before slow-roll conditions are no longer valid. Indeed, we 
start by considering the Lagrangian density in ( |l.lD and that both (j) and x will obey slow- 
roll conditions in this regime . In what follows and in conformity with the approximation 
scheme adopted, all quartic terms in the speed of and x will be neglected. 
The momenta conjugate to (j) and x ^-^^s, respectively 



dC 



7(/) and :- 



dC 

dx 



X 



Because the Lagrangian density is explicitly time- independent, it follows that the total 
energy is conserved 



E = j-) + \ + V{^) 



(A.l) 



Moreover, immediately before the beginning of the oscillations, we expect the time depen- 
dence of X to be parametrised by that of (p. Also, when the tachyonic instability sets in. 



the evolution of the fields will follow the linear trajectory (3^). In addition, writing for 
simplicity F((/>) = A + B(j)'^, with A = A{t) and B = B{t), following the evolution of the 
matter field, then from (|A.1|) 



E + 



f 



+ A + B<p'^ 



A B 

To proceed, we consider a "short" time evolution. At, so that — — > At and > At. 



A 

Redefining T :=£' + !// — A~ constant and keeping the leading order terms 



B 



P 



2BT(t? 



r(2 + /r) 
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Clearly, as we expect the inflaton field to be rolling down the potential, we take the negative 
root as the physical solution, so that 



1 



1/2 



BT(I? 

7^2 L 

P 



[T(2 + /T)]i/2 

for small 0. At this stage, we require T < and |/T| > 2, so that 

^ = C [1 + DcP^) , 

with 



(A.2) 



C 



D 



rn2 

\ T(2 + /T) 



B \T\ 

7^2 L 

P 



tan (c /D t + 13 /D 



This equation has solution 



where (3 is an integration constant which can be eliminated by requiring (p{t = tf) = 0, for 
some time tf, then 



tan 



C\fW {t-tf) 



(A.3) 



Manifestly, we are only able to study the evolution of (f) in such a regime in less than half a 

vr 

period, ie, for At < However, we can indeed check that 6^0, given the choice of 

2C^/IT 



the initial condition above. This kind of behaviour was interestingly obtained |3C] for the 
case of a finite throat, namely the Klebanov-Strassler throat, using the Hamilton- Jacobi 
method and neglecting any coupling of the inflaton with a matter field. 

Importantly, one would like to connect this behaviour with the oscillatory one when 
the bottom of the potential is finally reached. For small angles tan^ ~ sin 0/ cos ~ sin0, 
which is, apart from a shift in the phase, the same oscillatory behaviour for the background 
solution. The requirements also agree with the ones made earlier to get ( [A.3| ) and certainly 
restrict the energy values for which this study applies. 

B. Insight on Hill's equation 

Given Hill's parameters in to evaluate the solutions of Hill's equation of the form 

oo oo 



-2rzi 
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one usually applies Hill's method where recurrence formulae are arranged for the coefficients 

We will now use the well understood behaviour of unstable solutions of Mathieu's 
equation in order to infer what happens with respect to Hill's equation. We start by 
noting that in the edges of the resonance bands = 0; it turns out that for the most 
important instability band and for 62 , O4 < 1 



1 ± 



It 



6^2 



which, for 9\/62 <^ 1 (ie, well before the onset of reheating), essentially reduces to 



1 ± 



(B.l) 



From here, it also follows that, at lowest order, -/Oq' ~ I ± 62/2. Furthermore, from 
the instability parameter is given by 



sm- [ ifJ'k^J = A(0) sin ^6*0 - 



(B.2) 



where A(0) is the discriminant evaluated for fi = (following from Hill's method); A(0) 
is exclusively determined by the parameters ^2 and 64 and is approximately given by 



A(0) ~ 1 + cot (^/0^ ^ 



^ 

it can be shown that, approximately, this yields 

2 



+ 



1 - 00 4 - 00 



A(0) ~ 1 =F 



vr 



(B.3) 



Using ( B^ ) and (|B.3| ) , we find that from (lB.2| ) , we may compute an approximate expression 
for the instability parameter for Hill's equation as 



(B.4) 



This expression for is valid for 02,04 ^ 1- Also, we may approximate it even further so 
as to recover the expression for the Floquet exponent in the case of Mathieu's equation: 



1 



(B.5) 



This suggests that we evaluate the order of magnitude of ^2 and 04 in the early and 
late time regimes separately and then assess their relative importance in order to make 
further approximations and eventually convert Hill's equation into the standard Mathieu's 



equation. This motivates the study on section 3.2 
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